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. Abstract: We investigate whether interaction between massive neutrinos and 

, quintessence scalar field is the origin of the late time accelerated expansion of the universe. 

We present explicit formulas of the cosmological linear perturbation theory in the neutrinos 
probes of dark-energy model, and calculate cosmic microwave background anisotropies and 
qq ' matter power spectra. In these models, the evolution of the mass of neutrinos is deter- 

mined by the quintessence scalar field, which is responsible for a varying effective equation o 
of states; u e ff(z) goes down lesser than -1. We consider several types of scalar field po- 
^ . tential and put constraints on the coupling parameter between neutrinos and dark energy. 

By combining data from cosmic microwave background (CMB) experiments including the 
WMAP 3-year results, large scale structure with 2dFGRS data sets, we constrain the hy- 
pothesis of massive neutrinos in the mass- varying neutrino scenario. Assuming the flatness 
of the universe, the constraint we can derive from the current observation is ]T] m u < 0.45 
eV at la (0.87 eV at 2a) confidence level for the sum over three species of neutrinos. The 
dynamics of scalar field and the impact of scalar field perturbations on cosmic microwave 
background anisotropies are discussed. We also discuss on the instability issue and confirm 
that neutrinos are stable against the density fluctuation in our model. 
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1. Introduction 



After Type la Super Novae (SNIa) |Q] and Cosmic Microwave Background § observations 
in the last decade, the discovery of an accelerating expansion of the universe is a major 
challenge to particle physics and cosmology. Many models to explain such an accelerating 
expansion have been proposed so far, and they are mainly categorized into three: namely, a 
non-zero cosmological constant || , a dynamical cosmological constant (quintessence scalar 
field) [Q, 45], modifications of Einstein Theory of Gravity ||. One often call what drives 
the late-time cosmic acceleration as dark energy. 

While the existence of such dark energy component has become observationally evi- 
dent, the current observational data sets are consistent with all of the three possibilities 
above. The first observational goal to be achieved is therefore to know whether the dark 
energy is cosmological constant or dynamical component; in other words, the equation 
of state parameter of dark energy, w = P/p, is —1 or not. The scalar field model like 
quintessence is a simple model with time dependent w, which is generally larger than — 1. 
Because the different w leads to a different expansion history of the universe, the geomet- 
rical measurements of cosmic expansion through observations of SNIa, CMB, and Baryon 
Acoustic Oscillations (BAO) can give us tight constraints on w. Recent compilations of 
those data sets suggest that the dark energy is consistent with cosmological constant || . 

Further, if the dark energy is dynamical component like a scalar field, it should carry 
its density fluctuations. Thus, the probes of density fluctuations near the present epoch, 
such as cross correlation studies of the integrated Sachs- Wolfe effect |], || and the power 
of Large Scale Structure (LSS) |ji]], can also provide useful information to discriminate 
between cosmological constant and others. Yet, current observational data can give only 



poor constraints on the properties of dark energy fluctuations [|Tl], 12]. 

Another interesting way to study the scalar field dark energy models is to investigate 
the coupling between the dark energy and the other matter fields. In fact, a number of 
models which realize the interaction between dark energy and dark matter, or even visible 
matters, have been proposed so far [[b], 14, 15, [u], 17]. Observations of the effects of these 
interactions will offer an unique opportunity to detect a cosmological scalar field 13, 18]. 

An interesting model for the interacting dark-energy with massive neutrinos was pro- 
posed by Fardon, Kaplan, Nelson and Weiner [19, 20], in which neutrinos are an integral 
component of the dark energy. They consider the existence of new particles (accelerons 
with masses of about 10 -3 eV) that generate the dark energy, and the possible interactions 
between these particles and the known particles of the Standard Model. As far as quarks 
and charged leptons are concern, there are severe constraints on the strength of these in- 
teractions. However, neutrinos-weakly interacting very-low-mass particles- are exceptional, 
and cosmologically relevant couplings are still possible for them. This type has the scalar 
field adiabatically tracking the minimum of the effective potential, and, however, finally 
runs the adiabatic instability when the mass-varying neutrinos become non-relativistic, 
since the mass of scalar field must be much larger than the Hubble expansion rate for the 



adiabatic approximation to be valid [21 



The other type has the dark-energy provided by a quintessence scalar field evolving 



in a nearly flat potential, whose derivatives must satisfy the slow-roll condition. In our 
analysis, we adopt the latter type, in which the potential is somehow finely tuned so 
that the rate of scalar field evolution comes out to be comparable to the Hubble rate. 
Naturally couplings between the quintessence field to neutrinos leads to a back reaction on 
the potential, which was not taken into account in previous works [23, 24] with slow-rolling 
quintessence field, but we take into account fully the neutrino contribution to the effective 
potential in this work. We calculate explicitly Cosmic Microwave Background (CMB) 
radiation and Large Scale Structures (LSS) within cosmological perturbation theory. We 
first analyze constraints on a set of models, which has not been correctly done before, with 
cosmological observation data with WMAP3-year data and 2dFGRS data sets. 

In this paper, after reviewing shortly the main idea of the three possible candidates 
of dark energy and their cosmological phenomena in section II, we discuss the interact- 
ing dark energy model, paying particular attention to the interacting mechanism between 
dark-energy with a hot dark-matter (neutrinos) in section III. The evolution of the mass of 
neutrinos is determined by the quintessence scalar filed, which is responsible for a varying 
equation of states: uj e ff(z) goes down lesser than -1. Recently, perturbation equations for 
this class of models are nicely presented by Brookfield et al. p3j, (see also |24[]) which are 
necessary to compute CMB and LSS spectra. A main difference here from their works is 
that we correctly take into account the scattering term in the geodesic equation of neu- 
trinos, which was omitted there (see, however, the Erratum of Ref. ^5|). We will show 
and remark in section IV that this leads significant differences in the resultant spectra and 
hence the different observational constraints. In section V, we discuss three different types 
of quintessence potential, namely, an inverse power law potential, a supergravity poten- 
tial, and an exponential type potential. By computing CMB and LSS spectra with these 
quintessential potentials and comparing them to the latest observations, the constraints 
on the present mass of neutrinos and coupling parameters are derived. We also show that 
the equation of state w z= o can reach down to —1, which is consistent with experimental 
observations. In section VI, we discuss on the neutrino mass bound in our interacting 
dark-energy models. In appendix A, we show explicit derivation of geodesic equation in 
the Boltzmann equation which contains a new contribution from the time-dependent neu- 
trino mass and some formulas of the varying neutrino mass in early stage of universe in 
Appendix B. Finally the explicit calculation for the consistency check of our calculations 
in section III is shown in Appendix C. 



2. Three possible solutions for Accelerating Universe: 

Recent observations of SNIa and CMB radiation have provided strong evidence that we live 
now in an accelerating and almost flat universe. In general, one believes that the dominance 
of a dark-energy component with negative pressure in the present era is responsible for the 
universe's accelerated expansion. 

There are mainly three possible solutions to explain the accelerating universe. The 



Einstein Equation in General Relativity is given by the following form: 

G,u, = R^v-^Rg^v = 8-kGT^ + Ag^, (2.1) 

Here, G^ u term contains the information of geometrical structure, the energy-momentum 
tensor T^ u keeps the information of matter distributions, and the last term is so called 
the cosmological constant which contain the information of non-zero vacuum energy. By 
writing the Einstein equation with a flat Robertson- Walker metric, one can drive a simple 
relation: 

-0 + 3p) + -. (2.2) 



R 3 vr ' 3 

In order to get the accelerating expansion, either positive cosmological constant A (uj\ = 
P/p = — 1) or a new concept of dark-energy with the negative pressure (uu < —1/3) 
needs to be introduced. Another solution can be given by the modification of geometrical 
structure which can provide a repulsive source of gravitational force. In this case, the 
attractive gravitational force term is dominant in early stage of universe, however at later 
time near the present era, repulsive term become important at cosmological scales and 
drives universe to be expanded with an acceleration [26]. Also we can consider extra-energy 



density contributions from bulk space in so-called 'cosmological brane world' models, which 



can modify the Friedman equation as H oc p + p [27. pq [. 

In summary, we have three different solutions for the accelerating expansion of our 
universe as mentioned in the introduction. Probing for the origin of accelerating universe 
is the most important and challenged problem in high energy physics and cosmology at 
present. The detailed explanation and many references are in a nice review on dark energy 



1 29]. In this paper, we concentrate on the second class of solutions using the quintessence 
field. In the present epoch, the potential term becomes important than the kinetic term, 
which can easily explain the negative pressure with uj® ~ —1. However there are many 
different versions of quintessence field: K-essencc[]30], |3l|] , phantom j32|, quintom[^3|, and 
so on, and justifying the origin of dark-energy from experimental observations is really a 
difficult job. 

3. Cosmological perturbations in Interacting Dark-Energy with Neutri- 
nos: 

As explained in previous section, it is really difficult to probe the origin of dark-energy 
when the dark-energy doesn't interact with other matters at all. Here we investigate the 



cosmological implication of an idea of the dark-energy interacting with neutrinos [IS, 20]. 
For simplicity, we consider the case that dark-energy and neutrinos are coupled such that 
the mass of the neutrinos is a function of the scalar field which drives the late time accel- 



erated expansion of the universe. In previous works by Fardon et al. [19] and R. Peccei 



|2C], the kinetic energy term was ignored and potential term was treated as a dynamical 
cosmological constant, which can be applicable for the dynamics near present epoch. How- 
ever the kinetic contributions become important to describe cosmological perturbations in 
early stage of universe, which is fully considered in our analysis. 



3.1 Background Equations 

Equations for quintessence scalar field are given by 



+ 2^ + a 2 ^M = O, (3.1) 
d(p 

V eS (<f>) = TO) + TO) , (3.2) 

TO) = a" 4 / ^a/7+^TO)/(?) , (3-3) 

m v {$) = mp 1 ** , (3.4) 

where V(4>) is the potential of quintessence scalar field, Vi(4>) is additional potential due 
to the coupling to neutrino particles [jH], 34 1, and m v (<j)) is the mass of neutrino coupled 
to the scalar field, where we assume the exponential coupling with a coupling parameter 
(5. 7i is j j where the dot represents the derivative with respect to the conformal time t. 
Energy densities of mass varying neutrinos and quintessence scalar field are described 

as 

Pv = a" 4 J ^|^3 vV + a 2m lfo(q) , (3.5) 

3P U = a" 4 / rf3 g g 2 = f f q ) (3.6) 

7 (2tt) 3 vV + a 2 ™^ ^ v ; 

P* = ^ 2 + TO) , (3-7) 

p* = ^ 2 - to) ■ (3.8) 

From equations ( |3.5| ) and (3.6), the equation of motion for the background energy density 
of neutrinos is given by 

Q In hi, 

p v + m{ Pu + P„) = —Q^<t>{Pu - 3P V ) • (3.9) 

3.2 Perturbation equations: 
3.2.1 perturbations in the metric 

We work in the synchronous gauge and line element is 

ds 2 = a 2 (r) [-dr 2 + (<% + h i j)dx i daP] , (3.10) 
In this metric the Christoffel symbols which have non-zero values are 

r° 00 = \, (3-11) 

T% = ^S i:j + -hij + -kj , (3.12) 
nj = -J>) + \hij , (3-13) 
Tjfc = l^ihkaj + h aj> k - h jKa ) , (3.14) 



where dot denotes conformal time derivative. For CMB anisotropies we mainly consider the 
scalar type perturbations. We introduce two scalar fields, h(k,r) and r}(k,r), in k-space 
and write the scalar mode of h%j as a Fourier integral |35| 



hij(x, t) = Id ke 



where k = kk with k l h = 1. 



.3. jkx 



kikjh(k,T) + (kikj - r<5y)6r/(fe,r) 



(3.15) 



3.2.2 perturbations in quintessence 

The equation of quintessence scalar field is given by 

D(t>-V eS (4>)=Q 



(3.16) 



Let us write the scalar field as a sum of background value and perturbations around it, 
4>(x,t) = 4>(t) + 54>{x,t). The perturbation equation is then described as 



a a a Zi(x CLG) 



where 



dVi 



d 3 q de(q, 
(2vr) 3 dcp 



e{q,<f>) = V 'q 2 + a 2 ml(4>) 
de(q, 4>) a 2 m 2 (cf)) d In m v 



d(f) 



e(g,4>) d<p 



(3.17) 

(3.18) 
(3.19) 
(3.20) 



To describe 5 (^^j^J j we shall write the distribution function of neutrinos with background 
distribution and perturbation around it as 



f(x\ t, q, nj) = / (r, q)(l + V(x\ r, q, nj)) . 



Then we can write 



where 



dVi 



d 3 q d 2 e 
(2vr) 3 d<f> 2 



b + a 



d 3 q de 
(2vr) 3 ^ 



d 2 e a 2 ( dm u \ 2 a 2 m v ( d 2 m u \ a 2 m u ( de\ f dm v 



+ 



For numerical purpose it is useful to rewrite the equations ( 3.18| ) and ( 3.22| ) as 
dVi d In m v 



dVj\ _ d 2 h\m v 
~d$) = d(f> 2 



5<j>( Pu - 3P„) + 



dlnrriu 



(5p u - 35 P v ) . 



(3.21) 
(3.22) 

(3.23) 

(3.24) 
(3.25) 



Note that perturbation fluid variables in mass varying neutrinos are given by 



- 1 / d 3 q _ 4 f d?q de 



= a J (# e/o( ^ + a J WfT* Wo ' (3 - 26) 

MP* = «" 4 / /o(*)« - «" 4 / T^f 4SWo • (3-27) 



, (2tt) 3 e 7 (2vr) 3 e 2 

The energy momentum tensor of quintessence is given by 



T$ = HTM,* ~ \ (<T0,a + 2V {<!>)) % , (3.28) 



and its perturbation is 



2 

This gives perturbations of quintessence in fluid variables as 

1 dV 
~a~ 2 



' / " !r*6 a + 4>-> a bc\>^ + 2%-bc\> \ % . (3.29) 



-<5T ° = + ^-6<f> , (3.30) 



, ,~ 1 dV 



5P^ = -5T 73 = -^4>8<j) - —5<f> , (3.31) 



a 2 



k 2 ■ 

{p<t> + P<f>)0<p = ik'ST? = </>6<j> , (3.32) 
Sj = T] - <5jr^/3 = . (3.33) 
3.3 Boltzmann Equation for Mass Varying Neutrino 

One has to consider Boltzmann equation to solve the evolution of Mass Varying Neutrinos. 
A distribution function is written in terms of time (r), positions (x l ) and their conjugate 
momentum (Pi). The conjugate momentum is defined as spatial parts of the 4-momentum 
with lower indices, i.e., Pi = mUi, where Ui = dxij '(—ds 2 ) 1 ! 2 . We also introduce locally 
orthonormal coordinate = (t,r l ), and we write the energy and the momentum in this 
coordinate as (E,p l ), where E = \Jp 2 + m 2 . The relations of these variables in synchronous 
gauge are given by [35], 

P = -aE , (3.34) 
Pi = a(<% + hij)p> . (3.35) 
Next we define comoving energy and momentum (e, as 



e = aE = \J q 2 + a 2 m 2 , (3.36) 
qi = api . (3.37) 

Hereafter, we shall use (x l ,q,rij,T) as phase space variables, replacing f(x l ,Pj,T) by 
f(x l ,q,rij,T). Here we have splitted the comoving momentum qj into its magnitude and 
direction: qj = qrij, where n l rii = 1. The Boltzmann equation is 



in terms of these variables. From the time component of geodesic equation |3g] (see also 
Appendix A), 

~ = -T%P a P* ~ mg^m,, , (3.39) 



and the relation P° = a 2 e = a 2 \/q 2 + a 2 m^, we have 

dq 1 ■ „• „• 9 m dm dx' 1 , , . 

A = .. hijqn , n ,. a 2___. (3.40) 



Our analytic formulas in eqs.( [3.39 -3.40) are different from those of [^] and 24], since they 



have omitted the contribution of the varying neutrino mass term. We shall show later this 
term also give an important contribution in the first order perturbation of the Boltzmann 
equation. 

We will write down each term up to 0{h): 

d_l = dfo 5£ dfo 

or or or or 
dx 1 df _ q n i x j 
dr dx 1 e dx 1 ' 



dqdf _ f 2 m v dm v dx % 1 • { A dfo 
dr dq \ ° q dx 1 dr 2 U J dq 



drii df 



0(h 2 ) . (3.41) 



dr drii 

We note that §1 and J* are 0(h). 
3.3.1 Background equations 

From the equations above, the zeroth-order Boltzmann equation is 

dfo 



The Fermi-Dirac distribution 



* = AW = t5wr+I- <3 ' 43) 

can be a solution. Here g s is the number of spin degrees of freedom, hp and k-Q are the 
Planck and the Boltzmann constants. 

3.3.2 Perturbation Equations 

The first-order Boltzmann equation is 

^V*)"^£=£*-0. ,,44) 
e q 2 ocp dmq 



Following previous studies, we shall assume that the initial momentum dependence is 
axially symmetric so that depends on q = qft only through q and k ■ ft. With this 
assumption, we expand the perturbation of distribution function, ^f, in a Legendre series, 

ft, q, r) = ^H) £ (2^ + l)tt<(fc, q, r)P e (k • ft) . (3.45) 

Then we obtain the hierarchy for Mass Varying Neutrinos 

qi Q = — k^i + , (3.46) 

e 6 omq 

*i = i£fc(* -2* 2 ) + K , (3.47) 

*' = 7*(aTT* t --5^* t «) ■ (3 - 49) 

where 

K = _lg^^lnm^m/o 
3 e oq> o In g 

Here we used the recursion relation 

(£ + l)P m (/i) = (2* + l)/iP/M - • (3.51) 

We have to solve these equations with a g-grid for every wave number k. 

4. The impact of the scattering term 

In this section we mention the impact of the scattering term in the geodesic equation 
(Eq.( |3.39| )) due to mass variation on the CMB angular power spectrum. This term has 
been omitted in the recent literature [23, 24] (however, it was correctly included in the 



earlier work, see p6|). Because the term is proportional to ^ and first order quantity 



dx 

in perturbation, our results and those of earlier works [^, ^4|] remain the same in the 
background evolutions. However, as will be shown in the appendix, neglecting this term 
violates the energy momentum conservation law at linear level leading to the anomalously 
large ISW effect. Because the term becomes important when neutrinos become massive, 
the late time ISW is mainly affected through the interaction between dark energy and 
neutrinos. Consequently, the differences show up at large angular scales. In Fig. (|I]), the 
differences are shown with and without the scattering term. The early ISW can also be 
affected by this term to some extent in some massive neutrino models and the height of 
the first acoustic peak could be changed. However, the position of the peaks stays almost 
unchanged because the background expansion histories are the same. 

5. Quintessence potentials and Cosmological Constraints 



To determine the evolution of scalar field which couples to neutrinos, we should specify 
the potential of the scalar field. A variety of quintessence effective potentials can be found 




1 

Figure 1: Differences between the CMB power spectra with and without the scattering term in 
the geodesic equation of neutrinos with the same cosmological parameters. 



in the literature. In the present paper we examine three type of quintessential potentials. 
First we analyze what is a frequently invoked form for the effective potential of the tracker 
field, i.e., an inverse power law originally analyzed by Ratra and Peebles 



—^-\ (Model I) , (5.1) 

where Vq = M 4 = 3M^#^ stands for the vacuum energy, M p i is the plank mass and a 
is free parameter which will be constrained from observational data. 

We will also consider a modified form of V(<j>) as proposed by (3£| based on the condition 
that the quintessence fields be part of supergravity models. The potential now becomes 

Vtf) = M 4 e*PI*K (Model II) , (5.2) 

where the exponential correction becomes important near the present time as <p — > M p \. 
The fact that this potential has a minimum for <j) = \/ a /3M p \ changes the dynamics. It 
causes the present value of w to evolve to a cosmological constant much quicker than for 
the bare power-law potential |39|. 



We will also analyze another class of tracking potential, namely, the potential of ex- 
ponential type pQ| : 

V(<j>) = M 4 e~ a(l ^ ) (Model III) , (5.3) 

This type of potential can lead to accelerating expansion provided that a < In figure 
@, we present examples of evolution of energy densities with these three types of potentials 
with vanishing coupling strength to neutrinos. 




Figure 2: Examples of the evolution of energy density in quintessence and the background fields 
as indicated. Model parameters taken to plot this figure are a = 10, 10, 1 for model I, II, III, 
respectively. The other parameters for the dark energy are fixed so that the energy densities in 
three types of dark energy should be the same at present. 



5.1 Equation of State in Interacting Dark-Energy with Neutrinos 

In quintessence models, the scalar field 4> rolls down a self- interacting potential V(<p). The 
equation of state (EoS) is denned, 



_ 4?/2-v{4>) 

It has to satisfy the condition ljj, < —1/3 for the cosmic acceleration. Since lo^ > — 1 from 
eq.( |5.4| ), the oj^ < — 1 regime can not be realized by quintessence. However observations 
have shown exciting possibility of < —1, which causes the big-rip singularity problem 
by the phantom field at the end of the day. Various independent analysis of the Gold SNla 



data set [41] indicated it. Present situations of the observation tell us: First cosmic shear 
results from the Canada-France-Hawaii Telescope (CFHT) Wide Synoptic Survey provided 
a constraint on a constant equation of state for dark-energy, based on cosmic shear data 



alone: uj < -0.8 with 68% C.L. from the Deep Component of the CFHTLSpj. In the 



Supernova Legacy Survey (SNLS) [^] cosmological fits to the first year SNLS Hubble 
diagram gave u = —1.023 ± 0.090 ± 0.054 for a flat cosmology with constant equation of 
state when combined with the constraint from the recent Sloan Digital Sky Survey (SDSS) 



measurement of baryon acoustic oscillations. WMAP3 data [44] gave us two different 



results for different assumptions: when we assume flat universe including SNLS data, 

h0.07 
-0.09' 



uj = -0.97ln oL however if we drop prior of flat universe, WMAP + LSS + SNLS data 



provide ui = — 1.0Q2~^]^ and fifc = — 0.0241qqJ3. Still observations tell us the possibility 
of unexpected u < — 1. 
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Table 1: Equation of state in various quintessence potential models 
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Figure 3: Time evolution of the effective equation of state parameter w c s for models with potentials 
of inverse power law (Model I) and exponential types (Model III) . The parameters are fixed to the 
best fitting values except for those shown in the figure. The effective equation of state parameter 
can be smaller than — 1 at z > 0. 
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We summarize the possible values of EoS in various quintessence potential models in 
table In this section we show that < — 1 naturally arises when quintessence field 
interacts with neutrinos. Similar processes has been done before in the case of interacting 
dark-energy and dark-matter [16]. As pointed out by earlier works it is possible to have 
the observational equation of state w e g less than -1 in the neutrino-dark energy interacting 
models. The point is that any observer would unaware the dark energy interactions, and 
attribute any unusual evolution of neutrino energy density to that of dark energy. This is 
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seen as follows. Let us consider the recent epoch where the neutrinos have already become 
massive enough so that the energy density of neutrinos can be described as 

Pu = m v {cj))n v = m u ((j))n Ut0 /a 3 , (5.5) 

where n u ,0 is the number density of neutrinos at present time. One can decompose this 
into two parts as 

m "(0) iA „ fj. \_ u9 



Pu = m u ((j)o)n u o/a + — — - 1 m v (4>o)n v o/a , (5.6) 

\"V(0o) / 

and hence the Friedmann equation (neglecting baryon and photon contributions) 

H = — — (PCDM + Pu + P4>) 

8ttG ( a / m v {6) \ o \ 

= — ^— I (Pcdm.o + Pu,o) /a + I — -11 m v {(pQ)n vfi /a + p^j ■ (5.7) 

Therefore, while the first term in the above equation is regarded as a (total) matter density 
of our universe, the second and third terms comprise effective energy density which would 
be recognized as dark energy, 

/ m u {(f)) \ o 

/°eff = 7TT ~ 1 m A9o)nu,o/a° + p& . (5.8) 

\m u (4> ) J 

In observations one measures the equation of state of dark energy w e g defined by 

^ = -3H(l + w cS )p cS . (5.9) 
Here w e fi is related to the equation of state of quintessence through 

^eff = T~ ~~ , (5-10) 

1 — X 

m v {4>) A m^o)n„,o/o 3 

7-T- 1 i (5-11) 



which are derived from Eqs.(3.1), (3.7), and (|5.8| ). An example of time evolution of w e Q is 
depicted in Fig. (|3|). 

5.2 Time evolution of neutrino mass and energy density in scalar field 

For an illustration we also plot examples of evolution of energy densities for interacting case 
with inverse power law potential (Model I) in Fig. (Q). In interacting dark energy cases, 
the evolution of the scalar field is determined both by its own potential and interacting 
term from neutrinos. When neutrinos are highly relativistic, the interaction term can be 
expressed as 

( Pv - 3P V ) w -^{am v f p Vmzsslcss , (5.12) 

where Pi/ massloss denotes the energy density of neutrinos with no mass. The term roughly 
scales as oc a~ 2 , and therefore, it dominates deep in the radiation dominated era. However, 




scale factor 



Figure 4: Examples of the evolution of energy density in quintessence and the background fields 
in coupled cases with inverse power law potential (Model I). Model parameters taken to plot this 
figure are a = 1, j3 = 1, 3 as indicated. The other parameters for the dark energy are fixed so that 
the energy densities in three types of dark energy should be the same at present. 

because the motion of the scalar field driven by this interaction term is almost suppressed 
by the friction term, —2TC(p. The scalar field satisfies the slow roll condition similar to the 
inflation models, —27i.<j) ~ a 2 ^^(p u — 3P U ). Thus, the energy density in scalar field and 
the mass of neutrinos is frozen there. These behaviors are clearly seen in Figs. @) and 
. The derivation of the analytic expression for the behavior of varying neutrino mass in 
early time of the universe is given at Appendix B. 

5.3 Constrains on the MVN parameters and Neutrino Mass Bound 

5.3.1 Neutrino Mass Bounds from Beta Decays and Large Scale Structures 

In this section, we review in brief the present status of the determination of absolute 
neutrino mass from beta decay experiments and cosmological observation data. The exis- 
tence of the tiny neutrino masses qualifies as the first evidence of new physics beyond the 
Standard Model. The answer to the hot questions on (1) whether neutrinos are Dirac or 
Majorana fermions ? (2) what kind of mass hierarchy pattern they have ? (3) what are 
the absolute values of the neutrino masses, will provide us the additional knowledge about 
the precise nature of this new physics, and in turn about the nature of new forces beyond 
the Standard Model. There are three well known ways to get the direct information on 
the absolute mass of neutrinos by using: Tritium /3-decay experiment, neutrinoless double 
beta decay experiment, and astrophysical observations. 

The standard method for the measurement of the absolute value of the neutrino mass 
is based on the detailed investigation of the high-energy part of the /5-spectrum of the 
decay of tritium: 

3 H — ► 3 He + e~ + v e (5.13) 



1000 




1+z 

Figure 5: Examples of the time evolution of neutrino mass in power law potential models (Model 
I) with a — 1 and (3 — (black solid line), (3 = 1 (red dashed line), (3 = 2 (blue dash-dotted line), 
[3 = 3 (dash-dot-dotted line). The larger coupling parameter leads to the larger mass in the early 
universe. 

This decay has a small energy release (£"o — 18.6keV) and a convenient life time (T1/2 = 
12.3 years). Since the flavor eigenstates are different from mass eigenstates in neutrino 
sector, in general, electron neutrino can be expressed as 



5> 



(5.14) 



where v% is the field of neutrino with mass rrii, and U is the unitary mixing matrix. Ne- 
glecting the recoil of the final nucleus, the spectrum of the electrons is given: 



He 



dE' 



(5.15) 



and the resulting spectrum can be analyzed in term of a single mean-squared electron 
neutrino mass 



(mp) 2 =J2 m2 j\ U ej\ 2 = mi\U el \ 2 + m 2 2 \U e2 \ 2 + m 2 3 \U e3 \ 



(5.16) 
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If the neutrino mass spectrum is practically degenerate: mi — mi ~ 1113, the neutrino 
mass can be measured in these experiments. Present-day tritium experiments Mainz [56] 
and Troitsk|57] gave the following results: 



m\ = (-1.2 ± 2.2 ± 2.1) eV 2 (Mainz), 
= (-2.3 ± 2.5 ± 2.0) eV 2 (Troitsk). 

This value corresponds to the upper bound 

mi < 2.2eV (95%C.L.) 



(5.17) 
(5.18) 



(5.19) 
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Another useful method is by using the neutrinoless double beta decay. The search for 
neutrinoless double /3-decay 



(A,Z) — > (A,Z 2 ) +e~ + e 



(5.20) 



for some even-even nuclei is the most sensitive and direct way of investigating the nature 
of neutrinos with definite masses. In this process, total lepton number is violated and 
is allowed only if the massive neutrinos are Majorana particles. The rate of Qv(3(3 is 
approximately 

i = GWCQflB, Z) \M 0l/ \ 2 (m^) 2 , (5.21) 

J l/2 

where G 0l/ is the phase space factor for the emission of the two electrons, Mq u is nuclear 
matrix elements, and < man > is the effective Majorana mass of the electron neutrino: 



rrn 



(5.22) 



We can write eq.( |5.22| ), for normal and inverted hierarchy respectively, in terms of mixing 



angles and 



mn — m; 



7.9l|^10" 5 ey 2 , A a = ±(m 2 -m%) = ±(2.6 ±0.2) • 10" 3 eV* 



at the 3d level and CP phases as follows: 



(m ee ) 
{m ee } 



44mi + c 2 s 2 e^JA2 + m 2 + s 2 e^JAl + m 2 



s 2 mi + c 2 s 2 e^ JAl -Al + m 2 + c 2 s 2 e^ 



, (normal hierarchy)(5.23) 
(inverted hierarchy). (5. 24) 



However, Qv/3/3 decay have not yet been seen experimentally. The most stringent lower 
bounds for the life-time of O^/3/3-decay were obtained in the Heidelberg- Moscow^] and 



IGEX|U 76 Ge experiments: 



Ty 2 > 1.9 • 10 25 years (90%C.L.) Heidelberg - Moscow, 
Ty 2 ^ L57 • W 25 years (90%C.L.) IGEX. 



(5.25) 
(5.26) 



Taking into account different calculation of the nuclear matrix elements, from these results 
the following upper bounds were obtained for the effective Majorana mass: 



\mpp\ < (0.35 - 1.24) eV 



(5.27) 



Many new experiments (including CAMEO,CUORE,COBRA, EXO, GENIUS, MAJO- 
RANA, MOON and XMASS experiments) on the search for the neutrinoless double /3-decay 
are in preparation at present. In these experiments the sensitivities 



\m pp \ ~ (0.1 -0.015) eV 



(5.28) 



are expected to be achieved. More detail discussions will be appeared in the separated 
paper |60t|. 



Within the standard cosmological model, the relic abundance of neutrinos at present 
epoch was come out straightforwardly from the fact that they follow the Fermi-Dirac distri- 
bution after freeze out, and their temperature is related to the CMB radiation temperature 
Tcmb today by T v = (4/ ll) 1 / 3 T C mb with T C mb = 2.726 K, providing 

6C(3) _3 , , 

nv = TT^ T cmb, (5.29) 

where £(3) — 1.202, which gives n v ~ 112cm -3 for each family of neutrinos at present. 
By now the massive neutrinos become non-relativistic, and their contribution to the mass 
density (ft u ) of the universe can be expressed as 

n u h 2 = . (5.30) 

93.UeV v ; 

where S stands for the sum of the neutrino masses. In this relation, the effect of three 

neutrino oscillation is included [61]. We should notice that when obtaining the limit of 



neutrino masses one usually assumes: 

• the standard spatially flat ACDM model with adiabatic primordial perturbations, 

• they have no non-standard interactions, 

• neutrinos decoupled from the thermal background at the temperatures of order 1 
MeV. 

These simple conditions can be modified from several effects: due to a sizable neutrino- 
antineutrino asymmetry, due to additional light scalar field coupled with neutrinos p2[ . 
and due to the light sterile neutrino |33|]. However, analysis of WMAP and 2dFGRS data 



gave independent evidence for small lepton asymmetries p4| , 65], and such a scenario with 
a light scalar field is strongly disfavored by the current CMB power spectrum data ||66|| . 
We will not therefore take into account such non-standard couplings of neutrinos in the 
following. In addition, current cosmological observations are sensitive to neutrino masses 
0.1 eV < S < 2.0 eV. In this mass scale, the mass-square differences are small enough and 
all three active neutrinos are nearly degenerate in mass. Therefore we take the assumption 
of degenerate mass hierarchy. Even if we consider different mass hierarchy pattern, it will 



be very difficult to distinguish such hierarchy patterns from cosmological data alone [67]. 

After neutrinos decoupled from the thermal background, they stream freely and their 
density perturbations are damped on scale smaller than their free streaming scale. Con- 
sequently the perturbations of cold dark matter (CDM) and baryons grow more slowly 
because of the missing gravitational contribution from neutrinos. The free streaming scale 
of relativistic neutrinos grows with the Hubble horizon. When the neutrinos become non- 
relativistic, their free streaming scale shrinks, and they fall back into the potential wells. 
The neutrino density perturbation with scales larger than the free streaming scale resumes 
to trace those of the other species. Thus the free streaming effect suppresses the power 
spectrum on scales smaller than the horizon when the neutrinos become non-relativistic. 
The co-moving wave number corresponding to this scale is given by 

1 /2 

k nr = 0.026 nU 2 hMp C -\ (5.31) 



Cosmological Data Set 



CMB (WMAP-3 year alone) 
LSS[2dFGRS] 
CMB + LSS[2dFGRS] 

5) 

CMB + LSS + SNla 
CMB + LSS + SNla + BAO 

J) 

CMB + LSS + SNla + Ly-a 
CMB + LSS + SNla + BAO + Ly-a 



£ bound (2a) 



< 2.0 eV 

< 1.8 eV 

< 1.2 eV 

< 1.0 eV 

< 0.75 eV 

< 0.68 eV 

< 0.62 eV 

< 0.58 eV 

< 0.21 eV 

< 0.17 eV 



References 



Fukugita et al.[f71| 



Elgaroy et al.[72] 
Sanchez et al. 



Hannestad[74] 



Barger et al. 
Spergel et al. [ 76 1 



Goobar et al.[77] 

Seljak et al.|| 
Seljak et al.|| 



Table 2: Recent cosmological neutrino mass bounds (95% C.L.) 

for degenerated neutrinos, with almost same mass m u . The growth of Fourier modes with 
k > k nr will be suppressed because of neutrino free-streaming. The power spectrum of 
matter fluctuations can be written as 



P m (k,z)=P4k)T 2 (k,z), 



(5.32) 



where P*(k) is the primordial spectrum of matter fluctuations, to be a simple power law 
P*(k) = Ak n , where A is the amplitude and n is the spectral index. Here the transfer 
function T(k, z) represents the evolution of perturbation relative to the largest scale. If 
some fraction of the matter density (e.g., neutrinos or dark energy) is unable to cluster, the 
speed of growth of perturbation becomes slower. Because the contribution to the fraction 
of matter density from neutrinos is proportional to their masses (Eq. (p0|)), the lar Efcr 
mass leads to the smaller growth of perturbation. The suppression of the power spectrum 
on small scales is roughly proportional to f v |68[|: 



AP m (k) 

Pm(k) 



-8/,. 



(5.33) 



where f u = Q u /Qm is the fractional contribution of neutrinos to the total matter density. 
This result can be understood qualitatively from the fact that only a fraction (1 — f v ) of 
the matter can cluster when massive neutrinos are present p9| . 

Analyses of CMB data are not sensitive to neutrino masses if neutrinos behave as 
massless particles at the epoch of last scattering. According to the analytic consideration 
in |70|| ) since the redshift when neutrino becomes non-relativistic is given by 1 + z nr = 
6.24 -10 4 n v h 2 and 

z rec — 1088, neutrinos become non-relativistic before the last scattering 
when Q u h 2 > 0.017 (i.e. £ > 1.6eV). Therefore the dependence of the position of the first 
peak and the height of the first peak on £2 u h 2 has a turning point at £l u h 2 ~ 0.017. This 
value also affects CMB anisotropy via the modification of the integrated Sachs- Wolfe effect 
due to the massive neutrinos. However an important role of CMB data is to constrain other 
parameters that are degenerate with S. Also, since there is a range of scales common to 
the CMB and LSS experiments, CMB data provides an important constraint on the bias 
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parameters. We summarize some of the recent cosmological neutrino mass bounds in table 



5.3.2 Neutrino Mass Bound in Neutrino-Dark Energy Model 

As was shown in the previous sections, the coupling between cosmological neutrinos and 
dark energy quintessence could modify the CMB and matter power spectra significantly. 
It is therefore possible and also important to put constraints on coupling parameters from 
current observations. For this purpose, we use the WMAP3 |78|, fn| and 2dFGRS [g(J data 
sets. In figure |(| we demonstrate, as an example, the total neutrino mass contribution on 
the CMB angular power spectra with the inverse power law potential (Model 1). We can 
see that the deviations from observation data becomes severe when the total neutrino mass 
increases 



The flux power spectrum of the Lyman-a forest can be used to measure the matter power 
spectrum at small scales around z < 3 |81, 82]. It has been shown, however, that the 
resultant constraint on neutrino mass can vary significantly from ^ m u < 0.2eV to 0.4eV 
depending on the specific Lyman-a analysis used |83|| . The complication arises because the 
result suffers from the systematic uncertainty regarding to the model for the intergalactic 
physical effects, i.e., damping wings, ionizing radiation fluctuations, galactic winds, and so 
on p4|]. Therefore, we conservatively omit the Lyman-a forest data from our analysis. 

Because there are many other cosmological parameters than the MVN parameters, 
we follow the Markov Chain Monte Carlo(MCMC) global fit approach [85] to explore the 



likelihood space and marginalize over the nuisance parameters to obtain the constraint on 



parameter(s) we are interested in [86]. Our parameter space consists of 



P = (Q b h 2 , Q c h 2 , H, t, A s ,n s , mi, a, (3) , (5.34) 

where f^/i 2 and Q c h 2 are the baryon and CDM densities in units of critical density, H is 
the Hubble parameter, r is the optical depth of Compton scattering to the last scattering 
surface, A s and n s are the amplitude and spectral index of primordial density fluctuations, 
and (rrii,a,(3) are the parameters of MVN defined in section III. We have put priors on 
MVN parameters as a > 0, and [3 > for simplicity and saving the computational time. 

Our results are shown in Figs.(^) - (||). In these figures we do not observe the strong 
degeneracy between the introduced parameters. This is why one can put tight constraints 
on MVN parameters from observations. For both models we consider, larger a leads larger 
w at present. Therefore large a is not allowed due to the same reason that larger w is not 
allowed from the current observations. 

On the other hand, larger (5 will generally lead larger m v in the early universe. This 
means that the effect of neutrinos on the density fluctuation of matter becomes larger 
leading to the larger damping of the power at small scales. A complication arise because the 
mass of neutrinos at the transition from the ultra-relativistic regime to the non-relativistic 
one is not a monotonic function of j3 as shown in Fig.(||). Even so, the coupled neutrinos give 



10000 




Figure 6: The CMB angular power spectra for Model I with different neutrino masses: M„ = 0.3 
eV (upper panel) and M v — 1.0 eV (down panel). The solid line is the best fit for the model 
((a,(3) = (2.97,0.170)), the other lines are models with different parameter value of a and (3 as 
indicated. The points are WMAP three year data. 



larger decrement of small scale power, and therefore one can limit the coupling parameter 
from the large scale structure data. 

One may wonder why we can get such a tight constraint on /3, because it is naively 
expected that large (3 value should be allowed if fl u h 2 ~ 0. In fact, a goodness of fit 
is still satisfactory with large j3 value when Cl v h 2 ~ 0. However, the parameters which 
give us the best goodness of fit does not mean the most likely parameters in general. In 
our parametrization, the accepted total volume by MCMC in the parameter space where 
£l u h? ~ and (3 > 1 was small, meaning that the probability of such a parameter set is 
low. 
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Figure 7: Contours of constant relative probabilities in two dimensional parameter planes for 
inverse power law models. Lines correspond to 68% and 95.4% confidence limits. 



In tables ||-[|, we summerize our results of global analysis within la and 2a deviations for 
different types of quintessence potential with WMAP3 years data. We find no observational 
signature which favors the coupling between Mass Varying Neutrinos and quintessence 
scalar field, and obtain the upper limit on the coupling parameter as 

P< 0.46,0.47,0.58 {la); [1.11, 1.36, 1.53 (2a)], (5.35) 

and the present mass of neutrinos is also limited to 

fi A 2 oday < 0.0044, 0.0048, 0.0048 (la); [0.0095, 0.0090, 0.0084 (2a)], (5.36) 

for models I, II and III, respectively. When we apply the relation between the total 
sum of the neutrino masses M u and their contributions to the energy density of the 
universe: Vt u h 2 = M l/ /(93.1'leV), we obtain the constraint on the total neutrino mass: 
M v < 0.45 eV(68%C.L.) [0.87 eV(95%C.L.)] in the neutrino probe dark-energy model. 



The total neutrino mass contributions in the power spectrum is shown in Fig |10j, where we 
can see the significant deviation from observation data in the case of large neutrino masses. 
Before concluding the paper we should comment on the stability issue in the present 



models. As shown in |21, f22| , some class of models with mass varying neutrinos suffers 
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Figure 8: Same as Fig.(g), but for SUGRA type models. 

from the adiabatic instability at the first order perturbation level. This is caused by an 
additional force on neutrinos mediated by the quintessence scalar field and occurs when 
its effective mass is much larger than the Hubble horizon scale, where the effective mass is 
defined by = d 2 V e s / dcj) 2 . To remedy this situation one should consider an appropriate 
quintessential potential which has a mass comparable the horizon scale at present, and 
the models considered in this paper are the case p3|] . Interestingly, some authors have 
found that one can construct viable MVN models by choosing certain couplings and/or 
quintessential potentials [^7|, 88, 39 1. Some of these models even realises m e g 3> H. In 
Fig. (|IT|) , masses of the scalar field relative to the horizon scale m e fi/H are plotted. We find 
that m e s < H for almost all period and the models are stable. We also depict in Fig.(pT|) 
the sound speed of neutrinos defined by c;? = 5P u /8p u with a wave number k = 2.3 x 10 -3 
Mpc -1 . 
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6. Summary and conclusion: 

In summary, we investigated the dynamics of dark energy in mass-varying neutrinos. We 
showed and discussed many aspects of the interacting dark-energy with neutrinos sce- 
nario: (1) To explain the present cosmological observation data, we don't need to tune 
the coupling parameters between neutrinos and quintessence field, (2) Even with a in- 
verse power law potential or exponential type potential which seem to be ruled out from 
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Figure 9: Same as Fig.(M), but for exponential type models. 



Quantites 


Mean 


SDDEV 


lc-ranges 


2o"-ranges 


WMAP3 data (ACDM) 


a 


2.08 


1.35 


< 2.63 


< 4.38 




P 


0.38 


0.35 


< 0.46 


< 1.12 




n B h 2 [w 2 } 


2.21 


0.07 


2.15 - 2.28 


2.09 - 2.36 


2.23 ±0.07 


n CD Mh 2 [w 2 } 


11.09 


0.62 


10.52 - 11.68 


9.87-12.30 


12.8 ±0.8 


H 


65.97 


3.61 


62.30 - 69.37 


58.39 - 72.10 


72 ±8 


Z re 


10.87 


2.58 


9.81 - 12.15 


6.13-14.94 




n s 


0.95 


0.02 


0.94 - 0.97 


0.92 - 0.99 


0.958 ± 0.016 


A s [10 10 } 


20.66 


1.31 


19.38 - 21.92 


18.25 - 23.41 






68.54 


4.81 


64.02 - 72.94 


57.43 - 75.60 


71.6 ±5.5 


Age/Gyrs 


13.95 


0.20 


13.76 - 14.15 


13.59 - 14.40 


13.73 ±0.16 


n„ h 2 [w 2 ] 


0.36 


0.29 


< 0.44 


< 0.95 


< 1.93(95%C.L.)* 


T 


0.084 


0.029 


0.055 - 0.112 


0.031-0.143 


0.089 ± 0.030 



Table 3: Global analysis data within la and 2a deviations for the inverse power law type (Model-I) 
of the quintessence potential. (*)Using WMAP data alone, they found 1.93, but along with either 
SDSS or 2dFGRS galaxy redshift data, they have 1.83 or 0.97 at 95% confidence. 



the observation of lo value, we can receive that the apparent value of the equation of 
states can pushed down less than -1, (3) As a consequence of global fit, the cosmo- 
logical neutrino mass bound beyond ACDM model was first obtained with the value 
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Figure 10: Examples of the total mass contributions in the matter power spectrum in Model I 
(upper panel) and Model III (down panel). For both panels we plot the best fitting lines (green 
dashed), lines with larger neutrino masses M v — 0.3 eV (blue dotted) and M v = 1.0 eV (cyan 
dot-dashed) with the other parameters fixed to the best fitting values. Note that while lines with 
M„ — 0.3 eV can fit to the data well by arranging the other cosmological parameters, lines with 
M v = 1.0 eV can not. 

J2 m u< 0.45el/(68%C.L.) [0.87 eV{95%C.L.)]. 
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A. Boltzman Equations in Interacting Dark Energy-Neutrinos Scenario 

From the Lagrangian C = —m(4>)^ 

9tiv x 'fi%vi the Euler-Lagrange equation is given by 
d ( dC\ dC 



(A.l) 
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Figure 11: (upper panel): Typical evolution of the effective mass of the quintessence scalar field 
relative to the Hubble scale, for all models considered in this paper, (down panel): Typical evolution 
of the sound speed of neutrinos c 2 s — 8P u /8p u with the wave number k = 2.3 x 10~ 3 Mpc -1 , for 
models as indicated. The values stay positive stating from 1/3 (relativistic) and neutrinos are stable 
against the density fluctuation. 
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where 

dC x a 

— =P I1 = m(xn g m (A.2) 
dx^ y/-9iwXl*X v 

Therefore eq( |A.l| ) becomes 

1 d ( m {x») ^ , \ + m ^ 9ap,^ a x^ = _dm (A.4) 
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Quantites 


Mean 


SDDEV 


lc-ranges 


2o"-ranges 


WMAP3 data (ACDM) 


a 


6.19 


3.31 


0.10 - 7.77 


0.10 - 11.82 




P 


0.42 


0.42 


< 0.47 


< 1.36 




n B h 2 [w 2 ] 


2.22 


0.07 


2.16 - 2.29 


2.09 - 2.35 


2.23 ± 0.07 




11.10 


0.65 


10.47 - 11.75 


9.85 - 12.40 


12.8 ± 0.8 


#0 


65.37 


3.41 


61.84 - 68.73 


58.55 - 71.70 


72 ± 8 


Z re 


10.89 


2.62 


4.00 - 12.26 


4.00 - 14.78 




n s 


0.95 


0.02 


0.94 — 0.97 


0.92 — 0.98 


0.958 ± 0.016 


A S [W 10 } 


20.69 


1.32 


19.31 - 22.00 


18.20 - 23.32 




^ Q [10 2 ] 


67.90 


4.47 


63.74- 72.16 


57.59 - 75.02 


71.6 ±5.5 


Age/Gyrs 


13.97 


0.19 


13.78 - 14.16 


13.59 - 14.35 


13.73 ±0.16 


fi„ /i 2 [10 2 ] 


0.40 


0.32 


< 0.48 


< 0.91 


< 1.93(95%C.L.) 


T 


0.084 


0.029 


0.053 - 0.113 


0.028 - 0.139 


0.089 ± 0.030 



Table 4: Global analysis data within la and 2a deviations for the SUGRA type (Model-II) of the 
quintessence potential. 



Quantites 


Mean 


SDDEV 


lc-ranges 


2o"-ranges 


WMAP3 data (ACDM) 


a 


0.70 


0.42 


< 0.92 


< 1.41 




P 


0.50 


0.48 


< 0.58 


< 1.53 




n B h 2 [w 2 ] 


2.21 


0.07 


2.15 - 2.28 


2.08 - 2.34 


2.23 ± 0.07 


n CD Mh 2 [w 2 } 


11.10 


0.63 


10.48 - 11.72 


9.84-12.33 


12.8 ±0.8 


H 


65.61 


3.26 


62.37 - 68.70 


58.99 - 71.58 


72 ±8 


Z re 


11.07 


2.44 


10.07-12.35 


6.64 - 14.78 




n s 


0.95 


0.02 


0.94 - 0.97 


0.92 - 0.98 


0.958 ±0.016 


A S [W 10 ] 


20.73 


1.24 


19.48 - 21.95 


18.33 - 23.27 






68.22 


4.17 


64.38 - 72.08 


58.45 - 75.05 


71.6 ±5.5 


Age/Gyrs 


13.96 


0.19 


13.77-14.15 


13.61 - 14.36 


13.73 ±0.16 


Q„ h 2 [W 2 } 


0.38 


0.25 


< 0.48 


< 0.84 


< 1.93(95%C.L.) 


T 


0.086 


0.027 


0.058 - 0.113 


0.032 - 0.140 


0.089 ± 0.030 



Table 5: Global analysis data within la and 2a deviations for the exponential type (Model-Ill) of 
the quintessence potential. 



By using the relation ds = y 7 — g a px^x u d\, we obtain 



x^ . ...dx^ 

. = m(xH — — 

y/-g a f)&*xP ds 



and eq.(A.4) becomes 



m(x tl ) dx a dx@ 
~^~ 9a ^~ds~~ds~ 

dx? 



dm 

dx^ ' 

dm 
dx^ 



(A.5) 

(A.6) 
(A.7) 
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With simple calculation, finally we obtain the relations: 



-dT + T ^ P ~dT-~ 9 (A - 8) 

dP u 

'° + P a P P = -mg^rn, v . ( A.9) 



For /x = component, eq.( |A.9|) can be expressed as 



l|.(p0 ) 2 +r 0^ pap? = _ mg ^ m ^ (A.10) 



Since P° = g Po = a 2 e, each terms of the eq.( |A.10 ) are given by: 



-4n 2 „-4„^9 _-2rr„„2 , 2„_ ^ m 



First term = -2a~^Hq z + a" 4 ^- 1 - a~ z Hm 2 + a~~ 2 m— (A.ll) 

ar ar 



Second term = 2a~ i Hq 2 + oT 2 Plm 2 + a~ 4 -hi j q i q j (A.12) 



Third term = a 2 ra^^- (A. 13) 

ar 

Since the first term includes the total derivative w.r.t. comoving time, we obtain finally 
the eq. (ET4CD in Section III-C: 

dg 1 • , • 2 mdm dx % 

-r- = — haqn % n 3 - cf — ^— ——. (A.14) 

dr 2 JH q dx l dr K ' 

B. Varying Neutrino Mass in Early stage of Universe 

In the high-redshift region with z > 100, neutrinos behavior as like relativistic particles, 
even though they have non-zero mass. Since 1/e = [1 — a 2 m 2 /2q 2 ]/q we have 



(p-3i>„) = a~ 4 J 



d 3 q -' 2 — 



am 

Jo, 



(2vr) 3 e 

O -4 r f . , 2 /" ^ 



qdq(am) f - / —(am) / 



(2*) S 

u r " ■> . , / 4 4 



l-(amy-0(a*m*)} (B.l) 



(2vr) 2L 12 

For the massless neutrino case, the energy density is govern by 

/ \ -4 f d3( l f «" 4 7 4 Q \ 

(p,)ma SS fe ss = a J j^qf = j^y^ . (B.2) 

Then when we normalize it w.r.t. the energy density of massless neutrino: 

(p u -3P V ) 10 2 

— = — (am) (B.3) 

yyv /massless ' " 



Now we solve the equation of motion of quintessence field from eq. (|3,l| )with exponential 
type potential (Model III) for simplicity: 

7 . 2 dV 2 dV I 

= a 2 aV(0)-a 2 ^-(p u -3P u ) (B.4) 

In the relativistic neutrino case when they obey the slow-rolling condition: 

2H4> « -a JL (p u - 3 P v ) cc -p (B.5) 

Since the relative sign of (3 and <j> is always opposite, the slope of the varying neutrino 
mass become negative, this means that the primordial neutrino mass was decreasing when 
universe was expanded. 

C. Consistency check 

The form of k can be also obtained by demanding conservations of energy and momentum, 

i.e., demanding that V /t 5Ty +V ^ 5T„ = 0. Let us begin by considering the divergence 
of the perturbed stress-energy tensor for the scalar field, 



<« of- ■ 9 dV\ 9 / ■■ ■ , 9 d 2 V 

<T»= -a~ 2 cj) + 2H4> + a 2 — d u 5<p - a~ 2 5<t> + 2H6<fr + k 2 5<p + a 2 — 
\ d(p J \ dq> 

= ^W + ^W (c.i) 



where in the last line we used eqs. O and (p7|) . The diver gence of the perturbed stress- 



energy tensor for the neutrinos is given by, 



V M 6T{?= -bp - (p + P)diVi - m{5p + 5P) - h(p + P) (C.2) 
for the time component and 

V„ <5Tf= (p + P)vi + {p + P) Vl + 4H(p + P) Vi + d t P + djXi (C.3) 

for the spatial component. Let us check the energy flux conservation for example, starting 
with the energy flux in neutrinos (in fc-space): 

(pu + Pu)0„ = AxkaT* J q 2 dqqf (q)^i (C.4) 

where 9 U = ik l v u j. Differentiate with respect to r, we obtain, 

(p u + P v % + (p u + P V )0 V = 4irka- A J q 2 dqqfai>i - 4H{p v + P V )0 V (C.5) 



Let us consider the first term in the right hand side of the above equation. This gives 

= k 2 5P v - k 2 (p u + P v )a v + Unk 2 ^ 4 J q 2 



3 e 



2. Q 2 de 



+4irka 4 / q 2 dqqfoK 

where a is defined as (p + P)a = —(kikj — and expressed by the distribution 

function as 



8vr 

(Pi/ + Pv)&v = -yO 



' ' Q 2 dq?-fo(q)*2 



(C.6) 



Comparing eq. ( |C.5| ) with eq.( |C.3j ), we find that the divergence of the perturbed stress- 
energy tensor in spatial part for the neutrinos leads to 



9*V„ olf= -4vrA: 2 a- 4 



2 . <? 2 5e 
/; , — -vkk a I q dq-^ 



+ Airka 4 / q 2 dqqfoK 



■■> ... "•/-.hi i ........ I v /u" 

On the other hand, the divergence of the perturbed stress-energy tensor in spatial part for 



scalar field is, from eq.( C.l 



5 l V M <S7f = 



9 In m h 



(p v - 3P U ) = -kirk 2 5(t>a 



fq'dq^k- (C.8) 



These two equations imply that k shold take the form as eq. (|3.50 ). 

Next let us check the energy conservation. Density perturbation in neutrino is, (see 
eq.Q) 

d 3 q de 



&Pv = a 4 J ^-^36/0(9)^0 + 4 
By differenciate with respect to r, we obtain 

d 3 q 



bp v = —AH5p u + a 



(2vr) : 



(2vr) 3 d<p 
-4 / d3 Q 



(C.9) 



(2vr) : 



(2vr) 3 dr 



where 



/1 / 2 2 1 2 2 

e = [ria m + a m 



d In 



(2vr) 2 



0)A 



(CIO) 



d f de 



„ ,a 2 m 2 de „ , de d 2 e ■ 
-H — + 2H-- + —~d> 



(C.ll) 



Inserting eq.( |3.46| ) for \&o in the above equation, we obtain 

1 



-m{5 Pu + 5P V ) - {p u + p v ) 



+a 



d 3 q 



/( 



(2 7 r) 3 "' u 



d 2 



;5<t> + % 



de 



h(p u + Pv) 



(C.12) 



o 
o 

N) 
o 
o 
00 

o 
o 
o 
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Comparing with eq.( C.2 ), we find 



M 

V„ 57$= 



d 3 q , / <9 2 e _ , T <9e 

l/0 TTTo^ + ^O-^T 



(27r) 3 ' u V^ 2 



-i ; d 3 q de 



(C.13) 



which is found to be equal to — V M <JIg= —(5 



d(f> ' 
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